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(27 June 1996) 

In a recent Letter Jl| we presented the first numerical treatment of the full set of Boltz- 
mann equations for the evolution of an MeV Majorana r neutrino in the early Universe 
and concluded that mass limits from Big Bang nucleosynthesis were significantly weakened 
§ ' compared to previous investigations. 

Hi ! 

' An error in our numerical code unfortunately invalidates the results. As noticed by 
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several people [0 one should expect a strengthening rather than a weakening of earlier 
> ; m ass Ihnits based on the integrated Boltz-nann elation Bi because our .elusion of the 
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full neutrino spectra should reduce the annihilation of high-momentum neutrinos, thereby 
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permitting more r neutrinos to survive. This is indeed the case as shown in the revised Fig. 
1 m . Our neutrino distribution functions now differ less from kinetic equilibrium than found 
in [|l|]. A typical case for v T is shown in Fig. 2. Deviations by a factor of 2 appear frequently 
when compared to a kinetic equilibrium distribution, f v = [exp((E — fi)/T) + 1]~ , with 
the same number density and T = T 7 , as often assumed, whereas the distribution is close to 

X" 

that of a kinetic distribution with T and /x determined by the number- and energy densities. 

The resulting v e and distributions are still significantly heated relative to the standard 
case with a massless v T . For eV-mass v^, or v e this changes the present day contribution to 
the cosmic density to Q u h 2 = am^/93.03eV with a = 1 for a massless v T , and a =1.10(1.13), 
1.09(1.14), 1.03(1.05), 1.01(1.02) for u e (u^) for m Vr =5, 10, 15, 20 MeV (h is the Hubble- 
parameter in units of 100km s -1 Mpc -1 ). Later decay of v T can further increase the value 
of a. 

Fig. 3 illustrates the consequences for Big Bang nucleosynthesis in terms of the equivalent 
number of massless neutrinos, iV eq , needed to give a similar production of 4 He at a baryon- 
to-photon ratio rj = 3 x 10~ 10 . The revised results with and without inclusion of the change 



in the u e number density are in fine agreement with recent results based on the integrated 
Boltzmann equation by Fields, Kainulainen, and Olive ||. Including also the actual shape of 
the u e distribution to some extent compensates for the effect on the energy density, a result 
also found by Dolgov, Pastor, and Valle J7|||, though we disagree by several equivalent 
neutrino species with the total A^q found in || , and to some extent also with the differential 
changes in N eq quoted in [^]. 

In conclusion our revised results show that no MeV r-neutrino with mass below the 
experimental limit of 24 MeV is permitted unless more than 4 equivalent massless neutrino 
flavors become allowed by future observations of the primordial element abundances. 
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FIGURES 

FIG. 1. The relic number density of massive r neutrinos normalized to a massless species times 
mass, rm v . The solid curve is calculated using the full Boltzmann equation and all interactions. 
Dotted and dashed curves are adopted from || and 

FIG. 2. The distribution of r neutrinos of mass 15 MeV at a photon temperature of 0.92 MeV. 
The solid curve includes all possible interactions whereas the dashed-dotted curve only includes 
annihilations. The dotted curve is a kinetic equilibrium distribution with the same number density 
as that of the solid curve and T = T~, whereas the dashed curve has the same number- and energy 
densities as the solid, with T and fj, adjusted accordingly. 

FIG. 3. Equivalent number of massless neutrinos, shown as AN eq = N eq — 3. The dashed curve 
is without heating of u e , the dotted with only the electron neutrino number density changed, and 
the solid curve with the full distribution of v e . 
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